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We study two-dimensional U(iV) and SU(iV) gauge theories with a topological term on 
arbitrary surfaces. Starting from a lattice formulation we derive the continuum limit of 
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encoded in the theta term. In the topologically trivial cases the theta term gives only 
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gauge group on orientable surfaces it gives rise to a phase transition at 6 = ir, similar to 
the ones observed in other models. Using the equivalence of 2d QCD and a Id fermion 
gas on a circle we rewrite our result in the fermionic language and show that the theta 
term can be also interpreted as an external magnetic field imposed on the fermions. 
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1 Introduction 



Recently there has been a renewed interest in studying lattice models with a topolog- 
ical term. This was motivated mostly by the hope that the long-standing strong CP 
problem can be solved [lj. Recall that in four dimensional QCD one can in principle 
add a CP violating term, the so called theta term to the conventional Yang-Mills ac- 
tion. In nature however, strong interactions seem to respect the CP symmetry, so the 
coupling constant of this term (9) must be zero or at least very small. The question 
is why nature is forced to choose the special value 8 = 0. Up to now Monte Carlo 
studies of this model have been very limited mostly due to the following two difficul- 
ties. At first the theta term in Euclidean space-time is purely imaginary which makes 
the Boltzmann weight of some states negative. It means that the Boltzmann weight 
cannot be interpreted as a probability measure. This difficulty, however, has been 
largely overcome by some new simulation techniques developed in the last five years 
[0]. The other problem arising in Monte Carlo simulations is related to the fact that 
the theta term is of topological nature. In four dimensional QCD on the sphere the 
theta term is proportional to the instanton number, an integer that does not change 
upon continuous deformations of the field configuration. On the lattice however, the 
notion of continuity is completely lost, or in other words continuity does not impose 
any restriction on the deformation of the degrees of freedom at different locations. As 
a result there are no disconnected instanton sectors and the very definition of the topo- 
logical charge is ambiguous. A possible way of defining the topological charge is by 
picking a smooth "interpolating" field that is in some sense close to the given lattice 
configuration and assigning the topological charge of the interpolating field to the lat- 
tice configuration. This procedure is highly nontrivial and it is not clear what happens 
in the quantum continuum limit. An interesting question one can ask is whether the 
topological information encoded in the theta term is retained in the continuum limit. 

In this paper we shall study this question in two-dimensional U(iV) and SU(iV) 
Yang-Mills theories defined on surfaces of different topologies. In this case the topo- 
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logical term is proportional to the first Chern number of the underlying bundle. Our 
choice of the lattice theta term is the most naive, the only requirement it satisfies is 
that in the naive continuum limit it gives the Chern number. For U(iV) bundles on 
orientable surfaces the first Chern number can take any integer value depending on the 
"twist" in the bundle. On the other hand it is identically zero for the SU(iV) groups 
and also for U(iV) on non-orientable surfaces. Even in the latter cases, the same "theta- 
term" as for the nontrivial U(iV) case, can be included in the lattice action. With our 
choice of the lattice theta term it does not vanish for any of the above discussed cases, 
it does not necessarily give an integer topological charge on the lattice, but it has the 
correct continuum behaviour. We shall show that the topological information encoded 
in the theta term can be completely recovered in the continuum limit. This means that 
in those cases where the Chern number is trivial, the theta term completely decouples 
from the dynamics, giving only a constant shift of the free energy. On the other hand, 
for U(iV) gauge groups on orientable surfaces the theta term survives the continuum 
limit and it gives rise to a phase transition at 9 = tt, similar to the ones observed in 
the 2d U(l) gauge model [||] and in CP(iV) models at strong coupling [§J. 

The plan of the paper is as follows. In section 2 we calculate the partition function 
of 2d U(iV) and SU(iV) lattice Yang-Mills theories with a theta term using the standard 
group character ("strong coupling") expansion. In section 3 we derive the continuum 
limit of the action which will turn out to be the heat kernel with an additional contri- 
bution from the theta term. We interpret the continuum partition function in terms 
of the underlying topology and show the presence of a phase transition at 9 = it in 
the topologically nontrivial cases. Recently it has been shown that 2d QCD is equiv- 
alent to a one dimensional gas of fermions B. In section 4 we discuss our results in 
this language. In particular we show that the topological term can be interpreted as 
an external magnetic field imposed on this fermion gas. Finally in section 5 we sum- 
marise our results and discuss their possible relevance to more realistic models with a 
topological term. 
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2 The partition function on the lattice 



In this section we derive an expression for the partition function of 2d Yang-Mills 
theories with a theta term. We use the standard strong coupling character expansion. 
The expression of the partition function in terms of the character expansion coefficients 
is well known on arbitrary surfaces, see e.g. || [7j], so we shall just quote this result. 
We then calculate explicitly the character expansion coefficients for the specific lattice 
action that we pick. 

We use the conventional Wilson lattice gauge model; the degrees of freedom are 
gauge group elements associated to links of a two dimensional square lattice. The 
gauge group G can be thought of as either SU(iV) or U(iV) until some specific choice 
is stated. The action for a configuration \U{\ is given by 

s[u l ] = Y.s(u p ), (i) 

p 

where U p is the product of link variables around the plaquette p, s is a conjugacy 
class function on G (i.e. s{g^ 1 Ug = s(U),Wg, U G G) and the summation is for all the 
plaquettes. The partition function of a finite lattice consisting of A plaquettes is 

^) = £exp(-^y>(f/A (2) 

[Ui] \ 9 P J 

where g is the gauge coupling and the summation is for all the configurations. 

Since the Boltzmann factor on each plaquette is a conjugacy class function, it can 
be expanded in irreducible characters of G as 

exp = E d rMN/g 2 )Xr(U), (3) 

where 

Ar(iV/<? 2 ) = ^- J G dUexp (-^«(C0j X* r (U). (4) 

Here r labels the irreducible representations of G, d r is the dimension and \ r is the 
character of r. Substituting (§) into the expression of the partition function, the link 
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variables can be integrated out using the orthogonality of the characters. For a closed 
orient able surface of genus g we obtain |J 

Z A (g) = J2 d r- 29 ^(N/g 2 ) (5) 

r 

For non-orientable surfaces the summation is restricted to self-conjugate representa- 
tions and also there is an additional constant factor 0, which will be unimportant 
here, so we omit it. 

Up to this point we have not specified the form of the action, @ holds for any 
gauge invariant action. The information specific to the action is entirely contained in 
the A r 's. We now choose for the Yang-Mills part the usual Wilson action 

s w (U) = ±(trU + trlfl) (6) 

To make the analytic calculations feasible we choose 

Se {U)= l -{tiU-tiU ] ) (7) 

for the theta term. Obviously this have the desired continuum behaviour which can be 
easily seen in exactly the same way as the continuum limit of the Wilson term. While 
in the continuum limit s w (U) is proportional to the square of the field strength (up 
to a trivial additive constant), the term sg(U) goes linearly with the field strength. 
Notice also that the Euclidean theta term is purely imaginary as it should be. With 
these notations we parametrise the one plaquette action as 

s(U) = -coshe s w (U) + sinhe s e (£7), (8) 

where e is a real parameter that determines the relative weight of the two terms. This 
parametrisation of the action was chosen for convenience. If we wanted the relative 
weight of the theta term to be the larger one, the convenient choice would be to exhange 
the cosh and the sinh. The calculation would go in essentially the same way. 

In the reaminder of this section we explicitly calculate the integrals (f|) with this 
action in both the U(iV) and the SU(iV) case. 
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The integrand of (^J) is a conjugacy class function so we can rewrite it in terms of 
the eigenvalues of U, {exp(zai), exp(ia 2 ), exp(zajv)}- In terms of these variables the 
action is 



N 



s(a) = — y^(coshe cos a, — isinhe sin<x,- 



(9) 



The irreducible representations of U(N) can be labelled by N integers {n\ > ri2 >,••,> 
tin} = [n]. According to Weyl's character formula the character of this representation 
reads 

., e i<Xj(n k +N-k) 

X[n](a) = 3 , (10) 

where 

fill 



A (a) = dete ia ^ N ' k \ 

jk 



Finally the invariant measure on U(N) is given by 

N da 



dU oc J] ^ |A(3) 



m=l 



up to a factor independent of the en's. Putting this all together and introducing 

N 

f(a>) = — (coshe cos a — isinhe sin a) 
9 2 

as a shorthand notation we obtain 

d [n] A [n] (iV//) = const. J] T ^e^)A(a) (dete^+^r 



(12) 



(13) 



TV w 

const. Y\ / 
l=i J -* 



2tt 



jV 



(-1) CT JJ e^s^i) 

o-£5jv 3=1 



X 



A' 



e (-i)" n e 

r/gSjv m=l 



-ta m (n T) ( m )4-W-?7(m)) 



A? 



const, ^(-ir e n£^«xp 

/iGSjv crGSiv i=l 



f{ai)-iai(a(l)+n^ a{ i)-^a(T)) , (14) 



where the summations are on all the permutations of N elements. The terms of the 
last sum on a are easily seen to be independent of cr thus 

r7T da 



d [n] A [n] (N/g 2 ) = const. N\ det 

jk 



J(a)-ia(nj+k-j) 



7T 2vr 



(15) 



Now all that remains to be done is to calculate the Fourier expansion of e^ a ) 



n da 
— exp 



N 

— (sinhe cos a — i coshe sin a) — ima 

g 



da 

— exp 
-7T 2vr F 



N > 

— cosfa + it) 
T 



ima 



'I m (N/g 2 ). 



(16) 



Here I m is the modified Bessel function of order m. Substituting this into (|15D we 
obtain the final expression for the U(iV) character expansion coefficients 



A tT(») = T- d !f(« 



-e(nj-j+k) t 
J', 



nj —j+k 



(17) 



In the SU(N) case there are two differences compared to XJ(N). At first the ir- 
reducible representations of SU(N) are labelled by non-negative integers satisfying 
n± > ri2 >, ■■■ > Un = 0, which does not change anything in the previous calculation. 
Secondly for SU(iV) there is an additional Y^=\ a k = constraint on the eigenvalues 
that must be enforced by a delta function in the measure. If this delta function is 
written as 

N oo 

E e m ^ a \ (18) 



k=l 



the U(N) calculation goes through without any modification giving the final result 



A 



SU(N) 



(N/g' 



]T det L-^-j+k+m) In _ 



-j+k+m 



(N/g 2 ). 



(19) 



In the e — > limit the action is pure Wilson and the expressions we obtained give 
the correct character expansion of the Wilson action. Had we chosen to give a bigger 
relative weight to the theta term by exchanging the cosh and the sinh, we would have 
gotten the same result but with Bessel functions instead of modified Bessel functions. 



3 The continuum limit 

In this section we calculate the continuum limit of the character expansion coefficients 
to leading (nontrivial) order in g — > 0. In 2d gauge theories the convergence radius 
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of the strong coupling expansion is infinite and it makes sense to calculate the weak 
coupling limit of the terms in a strong coupling expansion. 

As g — > we want the coupling of the theta term to approach its continuum value 
9/2-k. Accordingly e has to be scaled in this limit as 

This can also be seen from dimensional analysis. For our purpose the only "interesting" 
part of the Aral's is the one that depends on the representation label. From now on we 
drop representation independent factors and whenever we do so the "~" sign will be 
used. Later on we shall comment on the information that might have been lost due to 
these omittions. 

Substituting the asymptotic representation of the modified Bessel functions 

e z _»£ 

I n {z) — ► 7 e 2z as z — > oo (21) 

\'2ttz 

and (p0| ) for e into (fL7|), for small g we obtain 

A W M) « -Ldete-^-J +k)2 ~^-J +k) } 
d[ n } i k 

1 / „2 N a \ 2 

w -- exp -A- Y nArii - 2j + 2 - -) det e-^^-^f*" 1 ) (22) 

In the last step the determinant was expanded as a sum over permutations, a factor, 
independent of the permutations was pulled out and the remainder was rewritten as 
a determinant. The last factor in fl22|) is a Van der Monde determinant which can be 
expressed as 

dete-^to-W* -1 ) = J| f e -£ (n >- j) - e -^ (nk - k A . (23) 
The last ingredient that we have not given yet is the dimension of the representation 

H 

d[n]~n(nk-k-nj+j). (24) 

j<k 

Making use of this and the asymptotic expansion 

f e -q*_ e -^ a + b 
in — ► Ini x + ... as x —>■ (25) 



one then obtains 

1 2 (a 2 N \ 

-j- det e-bto-M''- 1 ') ~ exp \-!L(N - 1) Vn,- (26) 

d [n] jk y N j=l J 

This completes the calculation of AFj in the continuum limit, giving 

/ 2 N n \ 

A^ N \g, 9) * exp f - A_ g^^. _ 2 j + N + 1 - -) J . (27) 

In the exponent for 9 = we recover the eigenvalue of the quadratic Casimir which 
gives the heat kernel action. 

The previous calculation goes through in the SU(iV) case without any modification, 
except that we have to do a little bit more work to evaluate the infinite sum in (|I"9|). 

oo .2 JV n 

Af n ] (7V W) « E exp-4£K. + m)(n, + m-2 J+ iV + l--) 

m=-oo /iV j=l 71 

oo 2 / V a \ 

= £ exp - A- iVm 2 + 2m£(n; - £-) . (28) 

m=-co ZiV y j=l Z7r J 

In the (7 — > limit this sum can be explicitly evaluated by approximating it with a 
Gaussian integral. After calculating this integral we obtain 

2" 



A su(N) ____ g 2 



exp 



2N 



^ , 6,1" N9 2 

i=i y=i 



GXp -2iV 



w fS^i^") 2 9 2 

X>,(^ - 2j + iV + 1) - V ^ Jy - iV^ 



(29) 



Again the # = case agrees with the character expansion of the heat kernel action. 

Notice that while in the U(iV) case 9 couples to the representations through their 
diagonal U(l) charge Y^=i n ji m the SU(iV) case 9 appears only in a representation 
independent term that can as well be dropped. This is exactly what we expected. Had 
we formulated the model in the continuum from the outset, the theta term would have 
been identically zero for SU(iV). The reason for this is that the continuum topological 
term is proportional to the integral of the trace of the field strength and the trace 
vanishes for SU(iV). 
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Let us now take a closer look at the U(iV) case. In the previous section we made a 
remark that on non-orientable surfaces only self conjugate representations contribute 
to the partition sum. In terms of the n/s the condition for a representation to be self 
conjugate is 

nj = -n N - j+1 Vj: l<j<N (30) 

It follows that for self conjugate representations X/JLi n j — and 9 cannot be seen in 
the continuum at all. Recall that in this case the topological charge must be identically 
zero because the second cohomology group of any compact closed non-orientable 2d 
manifold is trivial. Thus the integral of any closed 2-form and in particular that of 
the curvature must be zero. Indeed we do not expect the theta term to influence the 
continuum physics. 

Finally let us deal with the only case in which there is a nontrivial topological 
term in the continuum, i.e. the gauge group U(iV) on a compact closed orientable 
surface. In this case the configurations of the system are labelled by the irreducible 
representations of U(N), there is no additional constraint. The action corresponding 
to the representation [n] is given by 

1& E n >i -2j + N + l--) + (2g-2) \nd [nh (31) 

/iv j=i 7T 

where A is the number of plaquettes forming the surface. In the thermodynamic limit 
when A — > oo the last term becomes irrelevant and the partition sum is completely 
dominated by the ground state, the representation that minimises 

N Q 

8g [n] = E n A n i —2j + N + l ) (32) 

If 9 is varied, the ground state might change at some point (s) due to level crossing, 
which causes phase transitions. This phenomenon has already been observed in a 
variety of similar systems; 2d U (1) gauge theory with a theta term 0, U, the analogous 
Id spin model || and in 2d CP(N) models at strong coupling [|J. 

Let us now explore this possibility in the present model. For 9 = the action is 

N N N N 

soM = E njfa -2j + N + l) = J>? + Y, n j( N + 1 " 2 J) > E n ) (33) 
j=i j=i j=i j=i 
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The inequality can be easily verified by pairing the terms in the previous sum as 

N/2 

J2 [n k (N + 1 - 2k) - n N _ k+1 (N + 1 - 2k)] > (34) 

k=l 

and using the inequalities between the n^s. If N is odd then k goes only to and 
the remaining term vanishes. It is clear from ( |3"5| ) that at 9 = the ground state is 
given by the trivial representation n\ — ... — Un — 0. This is not surprising since it 
is well known that in the strong coupling expansion of the Wilson action the trivial 
representation gives the largest contribution for any value of the gauge coupling ||. 
Also from (^) we can see that 

se[n] > n Mi ~ -) = \ n i ~ + const - (35) 

which shows that by increasing 9 starting from zero we do not encounter any phase 
transition up to 9 = it. At this point however there is a phase transition; above % the 
representation rii = .... = n N = 1 becomes the ground state. Due to the periodicity in 
9 the same type of phase transition occurs whenever 9 becomes an odd integer times 
it. The free energy density around the 9 = tt phase transition is given by 

f(9) = < 9 < n 
f{9) = g 2 (l - ^ 9 > tt (36) 

We have to note that (|55|) gives only the part of the free energy that depends on the 
representation explicitly. The reason for this is that throughout the calculation we 
kept only the representation dependent terms of the free energy. Any representation 
independent but theta dependent neglected terms, however, should be analytic at the 
phase transition. Thus in the SU(iV) case, for example, we eventually would have 
needed to include an extra ^-dependent renormalisation counter-term in the free energy 
to make it ^-independent, as we expect it to be in the continuum. 
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4 The fermionic picture 

It is well known that 2d U(iV) QCD is esentially equivalent to a system of N fermions 
living on a circle || . It turns out that the result we obtained in the previous section has 
a nice interpretation in this language. The theta term produces an external magnetic 
field proportional to 9 and the phase transition we encountered at 9 = tt corresponds 
to an abrupt change of the ground state of the fermion system. This correspondance 
is possibly not very surprising since the theta term and the external magnetic field 
represent the same ambiguity in the quantisation of the two equivalent systems. This 
ambiguity is due to the non-simply connected nature of the classical configuration 
space. 

In this section we briefly summarise the fermion interpretation of 2d QCD from 
our point of view and discuss the relevance of the theta term. For a discussion of 
the fermionic picture starting from a continuum Hamiltonian formulation we refer the 
reader to [Q. 

In section 2 we obtained the partition function of 2d QCD using the Lagrangian 
description. Now we want to rewrite this in terms of the Hamiltonian in order to 
identify the (Euclidean) quantum mechanical system to which this partition function 
corresponds. In the Hamiltonian language the two coordinates on the surface corre- 
spond to time and one space direction. In order to simplify matters we first choose the 
simplest nontrivial surface, the cylinder. We take the compactified dimension to be 
"space" and time is measured along the axis of the cylinder. Using the gluing procedure 
described in section 2 and the form of the character expansion coefficients obtained in 
section 3 we can construct the partition function on a cylinder with arbitrarily chosen 
holonomies U\,Ui G U(N) on the boundary 

Z LT (9,Ui,U 2 ) = £e-^M xwCtfOxfniM, (37) 
N 

where L and T are the size of the cylinder in the space and time direction respectively. 
In the language of quantum mechanics the partition function of a system with fixed 
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configurations at t = and t = T is the (Euclidean) propagator 

(g 2 T| gi 0) = 5> 2 |j> e"^ T (38) 

where the |j)'s form a complete set of eigenstates of the Hamiltonian and the E/s are 



the corresponding energies. Comparing (|37|) and ( 38|) we conclude that in 2d QCD the 



energy eigenstates are labelled by the irreducible representations and the corresponding 
eigenfunctions are the characters. 

In this formulation the wave function is a function on the gauge group which might 
suggest that the classical configuration space is the group itself. This is however not 
quite right because the characters do not form a complete set in the space L 2 (U(N)), 
they span only the subspace of conjugation invariant functions on XJ(N). If we want 
to interpret the partition function of the cylinder as the propagator of a quantum 
mechanical system, the wave functions have to be restricted to be invariant under 
conjugation. This restriction looks more natural from the gauge theory point of view. 
At any fixed time t the only (independent) gauge invariant observable is the conjugacy 
class of the holonomy around the "space" circle at t. It follows that the classical 
configuration space is not the whole gauge group but only the set of its conjugacy 
classes. The most convenient parametrisation of this set is given by the phases of the 
eigenvalues of the group elements as described in section 2. In this way the classical 
configuration space in the U(iV) case can be identified with the iV-torus T N (a maximal 
torus in XJ(N)) factorised by the permutation group of N elements. This factorisation 
is essential because all the permutations of a given set of eigenvalues correspond to 
the same conjugacy class. The classical configuration space is thus equivalent to that 
of a system of N indistinguishable particles on a circle and the wave functions are 
symmetric functions on T N . 

There is however one more complication due to the nontrivial nature of the integra- 
tion measure defining the inner product on the Hilbert space. Recall that the invariant 
measure in terms of the above given coordinate system is given by flT^p. The inner 
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product is 

N da 



^^*(a) A* (a) A (a) #(a), (39) 



where A (a) is the determinant (11). It is now seen that by including the determinant 
A in the wave function the inner product is written in terms of the "usual" measure on 
T N . Since the original wave functions \P were symmetric the redefined wave functions 
^A will be antisymmetric. This shows that the particles we are dealing with are 
fermions. 

To see that it is actually the whole Hilbert space of N fermions that enters the 
expression for the propagator we have to look at the irreducible characters after the 
above redefinition 

X \n](a) A(a) = dete ia ^ nk+N ~ k) = dete ia ^ nk+ ^~ k) (40) 

jk jk 



and verify that they span the whole Hilbert space. (f5j) is essentially a Slater determi- 



nant composed of one-particle momentum eigenstates e iajPk with momenta 

N + l , . . 

Pk = n k A k. (41) 

For N odd it is easily checked that the irreducible representations are indeed in one- 
to-one correspondance with the states in the iV-fermion Fock space on a unit circle. 
If N is even the allowed momenta are half integers which give antiperiodic boundary 
conditions for the iV-fermion wave function. 



The wave function (4"0"|) describes an energy eigenstate with energy 



E [P\ = 7777 S *N 



2N 1 2N 



v / a \2 n , N+1 q x2' 



»( BY f 



J 2, 



(42) 



The first term is the energy of N noninteracting particles in the presence of a magnetic 
field with constant 8/2ir vector potential along the circle on which the particles live. 
The second term can be thought of as a trivial - though ^-dependent - shift of the 
energy that does not depend on the state of the system. For 6 < n the ground 
state is the usual fermionic ground state with all the states filled below the Fermi 



14 



surface. The Fermi surface in this case degenerates to the two points ±pf = ± 2^- ^ n 
summary, we have succeded in showing that 2d U(iV) QCD is equivalent to a system 
of A" noninteracting fermions moving on a unit circle. The correspondance between 
the amplitudes in the fermionic model and the gauge model is 

0T\aO) lcimion = A*0) A(a) (U0) T\U{a) 0} gauge (43) 

In this language the partition function on the torus can be obtained by taking 
periodic boundary condition in the time direction and integrating out the holonomy 
around the two identified boundary pieces 

Z to Ug) = Y,e- Tm (44) 

This is the partition function of the corresponding fermion gas at temperature 1/T. 
In the zero temperature limit only the ground state contributes to (f44|). At 9 = tt 



however the ground state becomes degenerate and above tt the new ground state will 
be the one with the momenta of all the particles shifted by one unit. This is the phase 
transition that we encountered in the previous section. 

Equation (K7I) gives the propagator in terms of momentum eigenstates. As a final 



check of our results it is instructive to rewrite the A^-fermion propagator using position 
eigenstates. This can be done by using the Poisson summation formula 

E /(P)= E / dxf(x)e 2 ~ (45) 



m=— oo 



The constraint p\ > p2--- on the momenta can be dropped at the expense of a factor 
AH and we also omit the representation-independent term in E[p\. Using the notation 
A [pi for the Slater determinant the fermion propagator reads 

(J3T\aO) e J — E, =1 (%-ft) = £ e -^E 3=1 fe-^) A[p] (a) 

Pj=— oo 



1 OO N r fOQ 



^ '"^ TLo 2 2 i • /n a , iTLa 2 a^ 

• P 2 lf-^+^i(27rm J -+a CT(j) -^ (j) + a ) 
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1_ (2Nn\> » e -fE^^V«rft) 



x x: (-ir«p 

cr,r;G5jv 



(46) 



This is the familiar form of the Euclidean propagator for a system of N fermions on a 
circle. Each term in the sum over m-configurations corresponds to the contribution of 
paths with a particular set of winding numbers [mi, mjv]- In complete agreement 

with our expectation the theta term gives an extra phase factor to each winding sector. 

Recall that if the configuration space of the classical system is not simlply connected 
quantisation becomes ambiguous. Paths belonging to different homotopy classes can 
get extra relative phase factors that belong to a representation of the first homotopy 
group of the classical configuration space. The ambiguity arises because there is no 
a priori criterium to select a particular representation JTIJ. The ^-dependent phase 
factor in (f46|) represents exactly this quantisation ambiguity, the representations of the 
homotopy group being labelled by 9. This concludes our discussion of the fermion gas 
interpretation of 2d U(iV) QCD. 



5 Conclusions 

In this paper we have shown a way of incorporating a theta term in two dimensional 
U(iV) and SU(iV) gauge theories. Starting from a lattice formulation we derived the 
continuum limit of the action resulting in a generalisation of the heat kernel action in 
the presence of a topological term. Our choice of the lattice theta term was the most 
naive, the only requirement it satisfied was having the topological term (first Chern 
number) as its naive continuum limit. In this way the lattice "topological charge" 
was not necessarily an integer and it was not necessarily zero in the topologically 
trivial cases. In spite of this, in the continuum limit we have been able to recover 
the topological properties of the system. In the topologically trivial cases when the 
topological charge was expected to vanish identically (SU(iV) gauge group on any 
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surface and U(iV) gauge group on non-orientable surfaces) the only effect of the theta 
term remaining in the continuum limit was a trivial constatant shift of the ground state 
energy. On the other hand in the topologically nontrivial case (U(iV) gauge group on 
orientable surfaces) the theta term remained to be coupled to the dynamics in the 
continuum, giving rise to phase transitions at 9 C = ±7r, ±37r, ±5tt... The same type 
of phase transition has been observed in CP(iV) models at strong coupling [|J. It is 
not very surprising that in the 2d gauge system even in the weak coupling limit we 
obtain an essentially strongly coupled system. This is due to the fact that in 2d QCD 
the strong coupling expansion has an infinite radius of convergence and we expect the 
strong coupling behaviour to persist down to any nonzero value of the coupling. 

Recently, based on Monte Carlo calculations, Schierholz argued that in the two 
dimensional CP(3) model and in four dimensional SU(3) QCD this phase transition 
also takes place but the critical value of 9 becomes small for weak coupling, eventually 
reaching zero in the continuum limit |l| . In this way 9 is forced to be zero if confinement 
is to be maintained and this gives a solution to the strong CP problem. From the 
argument given at the end of the previous paragraph it is clear that our results for 2d 
QCD give no direct hint concerning the phase transition in 2d CP(iV) or 4d QCD. Such 
a weak coupling transition must be directly related to the mechanisms responsible for 
maintaining confinement at arbitrarily weak coupling in four dimensions. We hope to 
address this question elsewhere. 

In the previous section we interpreted the continuum theta term in the fermionic 
formulation of 2d U(iV) QCD. In particular we have shown that in the fermion language 
the theta term corresponds to an external magnetic field imposed on the fermion gas. 
Very recently the connection between the gauge model an the fermion gas has been 
exploited to obtain some new results concerning the Calogero-Sutherland model [[l IJ . 
Our results might also find some application in this direction. 
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